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1. Introduction 

Developments in string theory suggest that matter and gauge interactions (described by 
open strings) may be localised on a brane, embedded into a higher dimensional space- 
time. Fields represented by closed strings, in particular gravity, propagate in the whole of 
spacetime. Such an idea provides an interesting framework to address theoretical issues of 
particle physics and cosmology, such as the hierarchy and cosmological constant problems, 
or the source of dark energy and dark matter. Several toy models were introduced (see for 
example ||l], ^ for early works) to address some of these problems. Some directly emerged 
from string theory setups [Q] while others were merely inspired by certain aspects of string 
theory. The possibility of large extra dimensions was considered in ref. |^] for an unwarped 
spacetime. It was later shown |^ that warping of five-dimensional spacetime could lead to 
localisation of gravity on the brane even though the extra dimension was of infinite proper 
length. Gravity in such models seems four-dimensional down to very small distances where 
higher (discrete or continuum) Kaluza-Klein modes take over. Later on several models 
were proposed where gravity was actually quasi-localised, i.e. four-dimensional at interme- 
diate scales and modified at very large distances (see for example [0, ^). Such models 
yield interesting possibilities for the origin of dark energy |jl^ , 'observed' in cosmological 
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experiments. However such models seem generically to suffer from ghosts [0] or strong 
coupling problems [12|. 

Let us step back to the case of localised gravity. As we mentioned above, gravity 
may be localised in the vicinity of a distributional, positive-tension, flat four-dimensional 
brane, embedded in a five-dimensional anti-de Sitter bulk, with an infinite extra dimension. 
The presence of a non trivial background is essential here. More precisely, the Einstein 
equations in this case admit a solution of the form 

ds^ = e2^(^) T]^^ dx^'dx'' + dz'^ . (1.1) 

The tension of the brane needs to be fine-tuned relative to the bulk cosmological constant 
in order to give four-dimensional Poincare invariance. This fine-tuning of parameters cor- 
responds to the usual cosmological constant problem in this context. Linear perturbations 
around this background show the existence of a single massless zero-mode, together with 
a continuum of massive Kaluza-Klein modes. The zero-mode is normalisable (because the 
volume element J e^'^dz is finite) and localised on the brane. This is interpreted as the 
usual four-dimensional graviton. Furthermore, the wave-functions of the massive modes 
are highly suppressed on the brane, so that their contribution to four-dimensional gravity 
at low energy may be negligible. This was further confirmed, at the linear level, by a 
rigorous derivation of Newton's law experienced by observers on the brane [p!^ ]. 

The Randall-Sundrum model Q and its basic features rely on assumptions, such as the 
localisation of matter on a distributional source (see however |^4|), a constant curvature 
bulk, and General Relativity to describe gravity in higher than 4 dimensions. One would 
like to know how generic these assumptions are. Since gravity plays such a key role in brane 
world models, it is certainly interesting to investigate how the usual features are modified 
by more general gravitational theories and backgrounds which are still compatible at some 
limit with Einstein's theory and fundamental principles, and are also in accord with string 
theory low energy effective actions. This is the direction we take in this paper. 

Toward this aim we add higher curvature contributions to the usual Einstein-Hilbert 
and cosmological constant terms in the bulk action. In four dimensions the Einstein-Hilbert 
term is the only term (apart from the cosmological constant) which yields at most second 
order derivatives in the field equations. In five dimensions this is no longer true and 
the quadratic Gauss-Bonnet term 

£gb = R^- ^RabR"' + R'^'^'Rabcd (1.2) 

also has this property. Any other quadratic term will produce higher order derivatives in 
the field equations, which will give rise to ghosts in the theory. Thus the Gauss-Bonnet 
combination is the only other sensible curvature term which can be included in the action 
for gravity (in five dimensions). 

The Gauss-Bonnet term is the Euler characteristic in four dimensions just as the Ricci 
scalar is the Euler characteristic of two dimensions. The similar properties of the two terms 
are direct consequences of this fact. In an arbitrary number of dimensions the most general 
action yielding second order field equations has been given by Lovelock [16| (see also [17| 
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for an elegant derivation using differential forms). Brane gravity with the Gauss-Bonnet 
term in a constant curvature bulk spacetime was first studied in |18|^. A careful treatment 



of the graviton fluctuations was given in ref. |19| (see also [^] for the full use of Lovelock 
theory in arbitrary dimensions). Cosmological evolution and consequences of such a setup 



have also been studied (see for example |2^, 22, 23, 24, |2^]). For some early applications 
to cosmology see | p7[ . 

Another generalisation is to add a scalar field, with a Liouville type potential in the 
bulk. This may represent the dilaton and/or moduli fields of compactified string theory. In 
this context, gravitational fluctuations around scalar brane backgrounds have been studied 
(see e.g. ^8|, and also for a smooth brane realisation). One of the phenomenological 
benefits of the inclusion of such a scalar field is an interesting reformulation of the cosmo- 



logical constant problem in brane world models [30|. One is then able to use the additional 
degrees of freedom to ensure the existence of a four-dimensional Poincare invariant solu- 
tion whatever the values of the brane tension and the bulk cosmological constant. Thus 
there is no fine tuning of the fundamental parameters of the theory. This is the so-called 
'self-tuning' mechanism. Unfortunately scalar field brane worlds with Liouville potentials 
generically suffer from naked singularities at finite proper distance from the brane, at least 
when sensible four-dimensional gravity is required [^]. This is true in particular for the 



self-tuning solutions [32|. The presence of the singularities re-introduces fine-tuning into 



the problem. This therefore represents a rephrasing, and not a solution, of the cosmological 



constant problem [33|. 



Close to these naked singularities, the original effective gravitational theory breaks 
down, because higher order terms in the Lagrangian cannot be neglected anymore. We 
are thus naturally led to consider higher order curvature terms and, for the same reasons, 
higher order scalar kinetic terms as well. In this paper we will consider terms which are 
up to quartic order in derivatives in the action. Such actions arise naturally in the context 



of low energy effective string theory |34|, and are known to lead to interesting properties 



in four-dimensional cosmology [^] as well as in black-hole physics |36]. Brane worlds with 



actions of this kind have been previously studied for a metrics of the form (LI) with A cx z 



The basic setup and particular solutions were analysed in detail in ref. [37|, and asymptotic 
behaviour of solutions and self-tuning was discussed in ref. |Q (see also |3^). In ref. pO|| , 
an exact solution was found with non trivial scalar field profile, no naked singularity and 
finite four-dimensional Planck mass. However, the 'self-tuning' mechanism was shown to 
fail since the solution was fine-tuned much in the same way as in the Randall-Sundrum 
model. In fact, as explained in section ^, in order to address this latter issue, general 
solutions have to be known analytically. The construction, analysis and stability study of 
such solutions is one of the aims of this work. 

The organisation of the paper is as follows: In the next section, we discuss the general 
action that we will consider. We demonstrate in particular that the origin of the scalar 
field (e.g. dilaton or moduli) fixes the coefficients of the higher order terms. We give the 
coefficients for toroidal Kaluza-Klein compactification. Field equations and boundary con- 



^see however |tl9|, [po| for the junction conditions 
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ditions for a conformally flat brane background ( |1 . 1[) are discussed in section ^. Next, 
in section we study the spin 2 linear gravitational perturbations around an arbitrary 



conformally flat background. We point out in all generality the essential differences of spin 
2 perturbations of Einstein-Gauss-Bonnet with respect to ordinary Einstein gravity. We 
address in particular the question of gravity localisation on the brane. It is also shown 
generically that stability requirements are far more stringent in higher curvature brane 
world setups. This essentially means that Einstein brane worlds can develop order a' 
instabilities in an effective action approach. Having discussed in a general context back- 
ground and perturbations, we move on to demonstrate how one finds solutions and their 
gravitational properties (singularities, localisation of gravity, stability etc). In section ^ we 
show how to obtain the general conformally flat solution to the field equations and their 
basic properties. In section |^ we apply our general analysis to the algebraically simplest 
case where the scalar field comes from the compactification of a flat sixth dimension. It 
is shown that the general maximally symmetric brane world solution is obtained from the 
corresponding six-dimensional black-hole solution. We concentrate on the cases without 
naked singularities in the bulk and study their stability and self-tuning requirements. It is 
intriguing that in the sole case where 'self-tuning' is actually possible (normalisable gravi- 
ton zero-mode without naked singularity in the bulk, nor brane tension tuning), a ghost 
appears in the bulk. The last section contains a summary of our main results and our 
conclusions. 

2. Action 

Consider the action describing a scalar field which is conformally coupled to gravity, written 
in the Jordan (or string) frame, with terms which are up to quartic order in the field 
derivatives 



This is the most general action of this order which obeys the rules of Lovelock gravity. 
This essentially means that the field equations obtained by variation of the scalar field (j) 
and the metric tensor g^^y are linear in second order derivatives, and do not feature higher 
order derivatives. This is an essential and natural hypothesis when studying backgrounds 
with boundaries and linear perturbations of them. In a string theory context the above 
frame is referred to as the string frame if (j) is the dilaton. Otherwise (for example if is a 
moduli field) it should be referred to as a Jordan frame. 

The constants a, b and c express the higher order terms for the scalar field and its 
interaction with the Einstein tensor. They will be given specific values later on, depending 
on the physical nature of the scalar field, uj is the Brans-Dicke parameter. The higher order 
terms couple via the dimensionfull constant a which is of dimension [length]'^. For string 
theory it is proportional to the string tension a' . Any terms other than those appearing 
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in the above action (|2.l|) yield higher than second order derivatives, which always result in 



the existence of ghosts |41], even around a flat background. 

Scalar-tensor theories of gravity are also studied in the Einstein frame, which is related 
to the Jordan frame by the conformal transformation gab — > e~'^'^^^gab- The Einstein frame 
action is 



M3 

Se = — 



j (fx^ |i? - ^ (3w + 4)5a09> - 2A 

£gb + 16a dactid^ct) - 166 G"'' dacpdbcp - 16c {da(l)d'' cp)^] } (2.2) 



where 

, 16-16 4, 8 , , 

a = a-3b + — , b = b-— , c = -2a + -b + c - — . (2.3) 

The identiflcations between the coefficients are easily obtained by comparing the field 
equations in the two frames (see Appendix). Written in this form the action more closely 
resembles the usual Einstein-Hilbert action. 

If the scalar field (p plays the role of the dilaton, the parameters a, b and c in eq. ( ^?l] ) 
are constrained in order to reproduce the scattering string amplitudes on shell [Q. For 
(J = — 1 this constraint reads, 

2a-26-c+l = 0. (2.4) 

The case a = b = c = 1 exhibits additional symmetry (higher-order extension of T- 
duality [^], as well as u-model conformal invariance p3[|). In the context of brane world 
models, it has been considered in ref. |38|, while the case a = 6 = 0, c = 1/27 has been 
discussed in detail in refs. |37, 

Another interesting case arises from toroidal Kaluza-Klein compactification. Start by 
considering the 5 -|- A^-dimensional action 

5^^+^) = ^ / d'^dX^y^^^^{R^'-'''^ - 2A + a£g+^)} (2.5) 

and take the extra dimensions to be flat and compact, with < < L and A,B = 
1 . . . N. Then for the metric ansatz 

dsl^^ = gabix)dx^dx' + e-4*(")/^r?^BrfX^dX« (2.6) 

we obtain the same field equations as eq. (|2.1D with 



uj = ^ ^ , a = {200 + 1)00, b = -uj, c = - {2oo + l)oo'^ . (2.7) 

This is the extension to higher curvature gravity theories of the usual Kaluza-Klein toroidal 
compactification. The scalar field in this case arises simply as the unique size or moduli 
of the A^ compact extra dimensions. For instance, in the string frame, A^ = 1 is simply 
oo = a = b = c = 0, and A^ = 2 is a; = —1/2, a = c = 0, 6 = 1/2. As A^ increases the 
extra terms in the action (^]^) are switched on, adding successive layers of complication. 



Interestingly the special case a = 6 = c = lof eq. (2.4) corresponds to A^ ^ oo and oo = —1 
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3. Field Equations 



Consider the conformally flat metric, 

ds"^ = e^^^'Uxf'dx^T^^^ + e^^^'Uz^ . (3.1) 

Four dimensional spacetime sections are Poincare invariant, and the single extra dimension 
(with coordinate label z) is warped. Note that we still have a gauge freedom relating the 
metric components A, B and (j)] for example B = corresponds to a Gaussian normal sys- 
tem. We shall use this gauge freedom to solve the field equations in convenient coordinate 
systems in sections ^ and ^. The full field equations are listed in the appendix. The (5,5) 
component of the generalised Einstein equation gives 

e^^A - 2uj(j)'^ + 6A'^ - 8(1)' A' + 

+4ae-2S (-3^'^ + 24A'^^' - SGbA'^cj)'^ + 4a {AA'^'^ + ^'^) - 6c(/.'^]) = (3.2) 

(note the absence of second derivatives). Furthermore a linear combination of the remaining 
two differential equations can be directly integrated to give 

g4A-2</-3B [e2B(2[a; + l]0' + A') 

-4a {3A'^ + 6A'2(/>' - 6b[2A'^^' + A'cj)'^] + Aa[2A'(l)'^ + (j)'^] - 4c(/>'^)] = -C , (3.3) 

where C is a constant of integration. It is interesting to note that for the Kaluza-Klein like 
choice of parameters (p. 7)), this equation factorises, 



g4A-20-B^2[a; + + A') (l - 4ae-^^[3A'^ + Gcj^V + 2a;(2a; + l)^'^]) = -C . (3.4) 

The integration constant C will acquire a physical meaning in the case of 1 dimensional 
toroidal compactification. The above system of equations (3.2-3^^) are independent and 
yield the full set of solutions of eq. ( p.l| ) for a conformally flat spacetime. 
The brane contribution to the action in the Jordan/string frame is 



J d^x^/^T{(t)) (3.5) 



where T is the brane tension. We can construct a brane world solution by joining together 
two spacetimes M.r and Ml- If we choose the coordinates so that z < in A^^, z > in 
Mr and z = on the brane, then the brane junction conditions are 



R 

= M-'^e^'t'T{(t)) (3.6) 

L 



[e-^{-4oj(l)' - 8A') + Aae-^^ {8A'^ - 24bA'^(l)' + 16a(^V'2 + (P'^S) - 8c(j)'^)]^ 



M- 



-3 20 



(3.7) 
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where [X]^ = Xn — Xi denotes the jump of X across the brane. It is interesting to note 
that a hnear combination of the above gives a similar equation to ds! 



[e-^(2(cj + l)(t)' + A') - 4Qe"2^ (3(^'3 + lA'^cj)') - 6b{2A'^^' + A'^'^) 

+Aa{2^r + 4>") - 4c<A'^)]f = (e'^ncP)) . (3.8) 



Thus if the integration constants in (3^) are such that [C]^ = (which corresponds to 
C = in the Z2-symmetric case), the RHS of ( |3.8D must vanish. This constrains the brane 
tension to be T{(j)) tx e~'^'^, which is the natural choice for the dilaton coupling to NS- 



branes in heterotic string theory. Such a case was studied in refs. |37, In ref. [40| in 
particular all requirements for 'self-tuning' (absence of naked singularities, finite volume 
element) were verified apart from 'self-tuning' itself. This was precisely due to the fact 
that the integration constant was set to be zero. Generically the [C]^ = solutions always 
require fine-tuning of the fundamental parameters. We will investigate in detail a case 
where C 7^ in section ^. 

4. Brane gravity and stability 

In order to investigate gravity as perceived by a four-dimensional observer and stability 
issues concerning the background solutions we will in this section consider perturbations 
around the general background ( |3.1| ). This will permit us to understand some important 
characteristics proper to the Einstein-Gauss-Bonnet theory that we are investigating. In 
this paper we will concentrate on the spin 2 perturbations which correspond to tensorial 
gravity on the four-dimensional brane world. We will see that the wave equation operator 
for the graviton modes in Lovelock gravity is closely analogous to that of Einstein gravity 



(see [51 1 for a recent discussion and references within). This follows from the fact that the 
two theories were constructed using similar requirements for their gravitational properties. 
The junction conditions will be however crucially different. 

Start by choosing an axial gauge for the perturbed metric (7a5 = 0): 

ds^ = e2^(^)(7?^^ + 7^^(x, z))dx^'dx'' + e^^^'Uz^ (4.1) 

and consider a perturbed scalar field (l){z) + f{x, z), with 7^jy and 99 small. The first order 
perturbation to the (//i/)-component of the gravitational field equation yields a rather 
complicated expression, 

d^{p{z)dz[j^u - Vtiul]) - w{z) (28" d(^,-i^)„ - U^^i^^ - d^d^^ - rif,y[d" d'^-iat, - □47]) 

= Kz)e'^ + 2r(0)) r,,,^ (4.2) 

where indices are raised and lowered with r/^,y, ^4 = 'q^^dfj^d^ is the flat four-dimensional 
Laplacian and /i, /2 and /s are some given background functions of z that we will not 
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need for what follows. What will be important however, are the functions w and p which 
are given by, 

w{z) = e2^-20+i? ^ Aaq{z)] + 4av'iz) , (4.3) 

p{z) = e4^-2<^-s (1 - 4aq{z)) , (4.4) 

where in turn 

q{z) = e-^^iA'^ - 4A'(I)' + 2^'^) , (4.5) 

We now project out the transverse traceless components of 7^,^ [Q (see also [^]). 
Define symbolically the operator -n^y = -q^y — dfj_dy/U\4^^ and consider 

1 



Itiv = \ ^^^,{a'^K)v - -^T^tiuT^aKj l"'^ ■ (4.7) 

This satisfies d^^p^y = and rf^^p^y = 0. Then, applying the projection ( p77| ) to the 
perturbation equation ([4.2|), we obtain the bulk perturbation equation for the tensor modes 

dMz)dz%u) + w{z)a4%y = . (4.8) 

Note that the form of the perturbation operator for ^f^y remains unchanged for a ^ 0. In 
particular it is free from derivatives of higher than second order. What changes however 
is the interaction with the background i.e. the background dependent functions w{z) and 
p{z) which are now of higher order in derivatives than when a = 0. This fact is spelt out 
by the definitions of q and v which are of first order in derivatives. Note in particular the 
presence of second order derivatives, v' , in the definition of w. This will yield important 
differences to the Einstein case for the perturbed junction conditions on the brane. 

Given the spacetime symmetries consider now a plane wave separation of variables: 

%yix,z) = c^y{x)i^m{z) With □4c;r,(x) = c^, (x) (4.9) 

where c^(x) is transverse-traceless and m? is the mass squared of the graviton mode as 
perceived by a four-dimensional observer. The scalar interaction wave equation ( |4.8D now 
reduces to 

{p{z)i,'^)' + rr?w{z)iP^=Q (4.10) 



valid everywhere in the bulk. ^Prom the explicit form of w and p in eqs. ( [4.3| ) and (4.4), 
the associated junction conditions at the brane are given by 

[e-^C]L = 4a[ge-^C + e-2^-^mVm(A' - 2<^')]f ■ (4.11) 

Furthermore, in order for the perturbed bulk to induce a well-defined metric on the brane, 
the modes ipm have to be continuous across the brane. 

In order to study the propagation of gravitons (existence of zero mode, stability, and 
orthogonality of eigenmodes) modulo the initial conditions supplied by the brane we now 
need to examine the boundary value problem as defined by ( 4.10| ) and ( 4.11 ). Such problems 



are generically treated in two ways; either as a Sturm-Liouville boundary problem, or by 
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coordinate transforming to a Schrodinger type equation. Let us follow these methods in 
turn pointing out the essential differences appearing here due to the higher curvature term 
corrections. 

If p and w are positive, the LHS of ( [4.10 ) is the usual 'self-adjoint' type differential 



operator with weight function w just like in standard (singular) Sturm-Liouville theory. 



The boundary conditions for (4.10) are however unusual and call for particular attention. 



For massless modes m? = 0, and Z2 symmetry, (4.11) tells us that = for z = just 



like for the a = standard Neumann boundary condition on the brane. For the 

massive modes however this is not the case, we have rather a mixed boundary condition 
and furthermore the mass of the mode in question appears in the boundary condition itself 
(4.11). This is precisely due to the presence of higher derivatives in w and is thus a feature 



of Gauss-Bonnet gravity^. We can however treat this in an analogous way by adding 
a suitable boundary term (total derivative) to the usual scalar product of the eigenvector 
functional space. Standard Sturm-Liouville considerations will then follow through. Indeed 
define 

{ct),i;)= (pijw dz + Aa['il;ct)v]f; (4.12) 

JBulk 

where not only w but also [v] has to satisfy positivity conditions. We have included a non- 
trivial boundary term which arises from the second order derivative appearing in w ([4.3|). 
Suitable Sturm-Liouville boundary conditions are imposed at asymptotic infinity. The 
problem ( p^HOiP is now self-adjoint The modes are then orthogonal, form a com- 



plete basis and the problem is free of tachyon modes. The particular boundary conditions 
imposed here ( [4.11 ) lead us to demand that [v] is also positive along with w and p. 



If on the other hand p is negative, the sign of the kinetic term in an effective action for 
■0m will also be negative. The solution will then have a graviton ghost in the bulk. If w and 
p have opposite signs the corresponding effective mass term will be negative. In this case 
the system may have tachyon modes and thus be classically unstable. Note that a negative 
value of [v] can also permit tachyon modes to exist, even ii w > (see section ^ for an 
example) . This is due to the mixed boundary conditions ( [4.11[ ) and essentially means that 
brane world backgrounds which are stable in Einstein gravity can develop an instability 
once we allow for a corrections. This amounts to a higher order curvature correction 
instability. This is the most important consequence of the mixed boundary conditions 
( [4.11 ) and is independent of the presence of the scalar field in the background. If p and 



w (and [v]) are negative the system will be classically stable although a graviton ghost is 
problematic at the quantum level. Note that the possibility of negative p or w (or [v]) only 
occurs when the a dependent terms are included in the theory. It should also be noted 
that even if there are no graviton ghosts or tachyons, there will be analogous issues with 
the scalar field perturbations (tachyons, ghosts and also radiative corrections), although 
we will not consider them in this paper. 

A necessary condition for gravity to be localised on the brane is that the zero mode 



^Not surprisingly similar boundary conditions also appear in the DGP model |^ with the inclusion of 
the induced gravity term on the brane. 
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solution of (|4.10| - ^n| ) be normalisable. This is equivalent to asking that 

w{z) \ ipo\'^ dz < oo . 



Bulk 



The zero mode solution of the bulk equation ( ^.lOP is 



I dz 

Vo = co + ci — , 

/Bulk P 



(4.13) 



(4.14) 



where cq and ci are constants fixed by the boundary conditions. For a Z2 symmetric brane 



world, ci = 0. If the zero mode ipo is constant the normalisation condition (4.13) amounts 
to having a finite four-dimensional volume element: 



Bulk 



w{z) dz , 



(4.15) 



which includes the a contributions. If the volume element ( [4.1 5|) is not finite and p > 
there will be no normalisable zero mode. However if we consider spacetimes where p has a 
different sign on each side of the brane, this argument may be bypassed. If w also changes 
sign the system may be classically stable, although it will have a ghost on one side of the 
brane. We describe an explicit example of this in section ^. 

It is also useful to recast the wave equation ( 4.10| ) in a Schrodinger form in the bulk 
{z / 0): 

d^ 



d^ + 



(4.16) 



with dx/dz = ^/w/p and = (wp)^^^ %l^rn- The effective potential is then given by, 



V{x) 



dx 



{wpyi^ 



Note that when a 7^ this equation is not valid on the brane since the coordinate x is ill 
defined there. The Schrodinger picture is however interesting when studying the large z 
behaviour. 



5. Conformally flat solutions 



We will now give a method for obtaining analytically the general solution of the bulk field 
equations (^I^^^l). If C = then e'^A' and e'^cj)' are constants and the differential 
equations reduce to algebraic equations, which are easily analysed. Such solutions have 



been investigated in [37, 38, 40 1. To deal with the more complicated case of C / 0, we 
begin by choosing our variable to be z = d(j)/dA (so in these coordinates, the brane is not 
necessarily at z = 0). 

The field equations now take the form, 



A + ^'2g-2B [Pq(^) _ 4ae-2^^'2Pi(z)] = 



(5.1) 
(5.2) 
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where 



Po = 6 - 8z - 2ujz^ , Pi = 3 - 24z + 366^^ - 16az^ - 2z^{2a - 3c) , 

Qo = l + 2{uj + l)z , Qi = 3 + 6(1 - 2b)z + 2z^{4a - 36) + iz^ia - c) . 
For A = 0, (5.1) is easily solved to give an expression for B 

B = lnA' + -ln^^ . (5.3) 
2 Pq 

Substituting this into ( |5.2| ) provides an expression for (j) — 2A in terms of z 

cj)-2A = ^ HPoQi - QoP,) + 1 In ^ - i ln(2CV^) . (5.4) 

Differentiating with respect to z and rearranging this expression gives 

A' = ^^^^_^^ ^{21n(PoQi -QoPi) + lnPo-31nPi} . (5.5) 

Integrating this gives A{z), which can then be substituted into eqs. (|5.3D and (5^) to 
find B{z) and (p{z). It is obvious that the choice of convenient a, b and c can simplify the 
relevant algebra. A similar approach can be used when A 7^ 0, although in this case the 
expression for B (5.3) will be a solution to the quadratic (5.1), which complicates all the 
subsequent steps. 

The properties of these solutions can be analysed by expressing the above polynomials 
in terms of their roots. Write Pq = 6(1 — uiz){l — U2z), Pi = 'iYli^ii^ — Viz), PqQi — 
QqPi = 15ni=i(l ~ Wiz). Some of the Ui,Vi,Wi may be zero, or complex. For example, if 
a = 6 = c = a; = 0all these parameters are zero except ui = 4/3, vi = 8, wi = —2. This 
case will be discussed in detail in section ||. 

Assuming that z = 2 is not a root of any of the three polynomials, the expressions 
(5^-5^) imply 

— — — In |1 — Uiz\ — - — — — In |1 — Vizl -\ — — — — In |1 — Wiz\ + est. 
l-2ui ' ' 4 ^ 1 - 2^;^ ' ' 2^ l-2wi ' ' 

i i i 

(5.6) 

lv-^1,, ,3v:-l,, 
(B = — > In 1 — Uizl > In 1 — Vjzl -\ — > In 1 — Wjzl + est. 

i i i 

(5.7) 

and ^ ^ 

P = In ^' - - ^ ln(l - Uiz) + ^Y1 ^^(^ ~ + ■ (^-^^ 

i i 

It is clear from the above expressions that the points z = 1/ui, l/vijl/wi (real values of z 
only) are singularities of some kind. We will refer to these as critical points. Note that the 
metric (and (j)) are only singular at the point z = ±00 if one of the critical points is zero. 

We will assume (for simplicity) that all the roots are distinct. In this case the behaviour 
of the metric and the curvature tensors near the singularities can be easily obtained from 



A 

4 
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the above expressions (^.3|-5.4). In the Jordan frame, we see that the points z = are 
curvature singularities which are a finite proper distance away from other nearby points 
(i.e. / e^dz does not diverge as z ^ l/^^i). On the other hand, the points z = l/uj, are 
aU at infinite proper distance from other points, and are not curvature singularities. They 
therefore correspond to spatial infinities. Hence all the roots of Pqi P\ and P^Qx — QoPi 
are either naked curvature singularities or spatial infinities. A slightly different analysis 
is required for 1/z 0, but the same results apply. The above expressions actually give 
several solutions; one for each range of z which does not contain any coordinate singularities, 
and for which B ( |5.8D is real. 

The situation is similar in the Einstein frame. The critical points are curvature sin- 
gularities at finite proper distance if the corresponding critical point has 1/2 < Ui < 2/3, 
Vi < 1/2, Vi > 1 or Wi > 1/2. Otherwise they are spatial infinities. The difference between 
the two frames arises because 4> — > ±00 at the critical points. 

We are interested in regions of spacetime which can be used to construct brane worlds 
with localised gravity, but which do not have any problems with naked singularities. If the 
zero- mode ipQ is constant as usual, the condition ( 4.131 ) for gravity to be localised on the 
brane is equivalent to require that the volume element (4.15) is finite. By substituting the 
solution (5.6-5!8|) into eq. ( 4.15 ), we can determine whether it converges near each of the 
critical points. If it does then we can construct a suitable brane world using such a region 
of spacetime. The volume element ( 4.15 ) will converge near a critical point if it satisfies 
1/2 < Wi < 1, Vi > 2, Vi < 1/2 or 1/2 < Ui < 2/3 as appropriate. Thus in the Jordan frame 
there are ranges of Ui and Wi for which regions of spacetime without curvature singularities 
and finite volume element can be found. However it should be noted that if Ui,Wi > 1/2 
(or Vi < 1/2) then the string coupling, e"^, diverges at the corresponding critical point. 
Furthermore the time taken (as perceived by an observer on the brane) for a photon to 
reach the brane from this point (J e^~^dz), is finite if Wi > 1/2, Wi < 0, Vi > 2, Vi < 1/2 
or 1/2 < lij < 2/3. Hence the region of strong coupling will be visible from the brane. In 
the Einstein frame the conditions for a finite volume element are the same as those for a 
naked singularity. 

If A 7^ a similar approach can be used. Asymptotic behaviour of the solutions can be 
obtained near the critical points, and this can be used to determine if the volume element 
is finite or if there are naked singularities. As for A = there are no solutions with 
finite volume element and no curvature or string coupling singularities. Similar ideas have 
previously been used to analyse the a = b = c = —uj = 1 case . 

Hence, for C 7^ 0, it does not generally seem to be possible to obtain brane models with 
finite volume element and without naked singularities. This is not true for C = |38, 40 1, 
but the brane tension is then fine-tuned, as explained in section ^. 

There are several cases for which the above analysis does not apply. If the polynomial 
q ( [4 .51 ) which appears in the definition of the volume element ( |4.3| ) shares a root with Pi, 
the volume element will converge for a wider range of parameters. Similarly if 1 -|- 4ag 
shares a root with one of the other polynomials. The situation will also be changed if any 
of Pq, Pi, or PqQi — QqPi share a root. However these situations will not generally arise 
without unnatural fine-tuning of the parameters oj, a, b and c. An exception is the point 
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z = oo. For many of the Kaluza-Klein like choices of parameters ( |2.7D this critical point 
corresponds to a root of more than one of the polynomials. However each of these cases 
has similar properties to the non-degenerate critical points discussed above. 



6. Six-dimensional Kaluza-Klein case 

Let us now focus on the case where 4> is a moduli field of a compact sixth dimension. In 
the Jordan frame, putting = 1 in the action ( |2.5D corresponds toui = a = b = c = 0. 
This is obviously the simplest case to consider for the action We will now construct 

the general solution for these parameters. 

6.1 General bulk solution 

We start by considering the 6-dimensional Gauss-Bonnet anti-de Sitter black hole solu- 



tion |4g]. The solution reads, 



dsl = -V{r)dT^ + r^dxl + (6.1) 



where V{r) = k + r^{l ± U)/{l2a) with 



Here, dx| is an Euclidean space of constant curvature k, and Ac = —5/{12a). Note in 
passing that if A = Ac the field equations have a degenerate solution with A' = 1/ ^/\2a 
and (f) completely arbitrary. This corresponds to the Class I solutions found in ref. pi] ]. 

The integration constant /i is related to the 'AD energy' of the solution. Unlike ordinary 
Einstein gravity (a = 0) there are two branches of solutions and both vacua are classically 
stable [^] to small perturbations. The solution corresponding to the (— ) choice of sign 
has a well defined a ^ limit, whereas the second solution exists only if a 7^ 0. For lack 
of a better name we will call the (— ) branch the Einstein branch, and the (-I-) branch the 
Gauss-Bonnet branch. Note that particular attention has to be given to the 'AD energy' 



of the Gauss-Bonnet branch as was shown in ref. (see also |48] for a Hamiltonian 
approach). We emphasise that the Gauss-Bonnet branch does not have an equivalent 
solution in Einstein theory. 

Let us focus on the planar case k = 0. For both branches there is the standard 
curvature singularity at r = 0. If (1 — A/Ac)/i < then there will be an additional 
singularity at r = r,,, = (A/Ac — 1)^"^^^/U^/^, where U = Only the solution with the lower 
choice of sign has an event horizon, and then only if A/i < 0. The horizon is located at 
r = rh = (Ac/A)^/^^^/^. Otherwise there is no black hole but merely a naked singularity. 

A conformally flat, five-dimensional scalar field solution is obtained from the black 
hole (|6.lD by the analytic continuation t = ix^ and Xi oc iT, comparison with ( |2.6D , and 
then compactification of the sixth dimension. The scalar field solution in brane background 



( |3.1| ) reads, 

A = In - + A , (6.3) 

n 
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1 r^(l±U) 

B = -^lnr^{l±U) + ^lnl2a , (6.5) 

where the parameters Af,, (ph and are constants which wiU be convenient when discussing 
boundaries, and 




Ub = U{n) = 

Note that the parameter n in the above solution is proportional to the constant C in 
Taking Af, = {2/3)4>h the metric in the Einstein frame reads, 

If the original vacuum solution has an event horizon, it will be transformed into a cur- 
vature singularity for the scalar field system. This is a common phenomenon of conformally 
flat scalar field spacetimes obtained in this way. Of course a naked singularity signals the 
endpoint of the validity of a particular solution. There are a total of 7 distinct types of 



spacetimes which can be obtained from the metric (6.1) in this way. They can be classified 
according to the choice of sign in the definition of V{r), the sign of ^u, and the value of A. 
Most of them have naked curvature singularities which restrict the allowed range of r. 

For the Gauss-Bonnet branch, if A < Ac (in which case must be positive), the metric 
is only well-defined in the coordinate range < r < r*. This solution has two naked 
singularities which are a finite proper distance apart. If A > Ac the solution is well-defined 
for r > r=K if /u < 0, and for r > if ^ > 0. These solutions have one naked singularity and 
are of infinite proper distance i.e. infinite in the extra dimension. 

For the Einstein branch we have more possibilities due to the additional singular point 
r = rh- The metric is only well-defined for r/^ < r < r,,, if A < Ac (and /i > 0), and 
for r=K < r < r/i if A > (in which case fj, must be negative). Again we have two naked 
singularities separated by finite proper distance. When Ac < A < the solution is well 
defined for r > if fi > and for r > r* if /i < (or just simply r > if ^ = 0). If A = 
for the Einstein branch, must be negative and the metric may be defined for r > r^,, but 
has a naked singularity at infinity (r oo). 

Finally for A = Ac, there is only one branch ^, 1 ±11 being replaced by 1 -|- fl/r^^'^. 
The solution is well-defined for r > (—p,)"^^^ if /i < and for r > if /i > 0. 

For comparison with the previous section, the three A = solutions: Gauss-Bonnet 
branch with < 0, Gauss-Bonnet branch with fi > 0, and Einstein branch with /i < 
respectively correspond to the three ranges — oo < z < —1/2, —1/2 < z < 1/8 and 
3/4 < z < cxD of the coordinate z = d(j)/dA. The points z = —1/2,3/4 are spatial infinities 
(r = cxd), while z = 0, 1/8 are the curvature singularities {r = r^:,0 respectively). 

An infinite (proper distance) brane world solution can be constructed by taking two 
halves of infinite bulk spacetimes (not necessarily identical) and 'gluing' them together. The 

^In D = 5 this case is related to Chern-Simons gravity (see for instance |M9|) 
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brane acts as a boundary for the bulk spacetime on both sides permitting the ehmination of 
one curvature singularity. We see immediately that spacetimes with two naked singularities 
(which are also of finite proper distance) are unsuitable for constructing singularity free 
brane worlds. Thus we can construct singularity free brane world models using the any of 
the A > Ac Gauss-Bonnet branch solutions or the Ac < A < Einstein branch solutions, 
provided that we keep the appropriate singularity free half of the bulk spacetime. Note 
that we keep asymptotic infinity r — > oo which represents proper infinity too. We take the 
brane to be positioned at r = r;,. The constants A^, and (ph have to be the same on both 
sides of the brane, in order to guarantee a well-defined induced metric as well as continuity 
of the scalar field. In contrast, the parameter /i, or equivalently Ub, may take different 
values on each side. 

Before moving on to boundary conditions we note in passing that if the black hole 



horizon is not flat then analytic continuation (see for instance |50] for an a = example) 
gives in the Einstein frame, 

with (f> = —^lnV{r) + Notice then how the curvature of the horizon after analytic 
continuation translates into the de Sitter or anti-de Sitter expansion of the brane world 
sections. The metric ( |6.7D represents the bulk background for a brane world of constant 
curvature. Having obtained the full set of spacetime solutions we now can move on to the 
boundary conditions. 

6.2 Brane tension and parameter tuning 

For a Z2-symmetric brane world, the solution ( |6.3| -|6.5|) is defined up to 2 integration 
constants, namely (j)b and fj, (the constant Ah is not physical and may be taken to vanish), 
and 2 possible branches. Applying the junction conditions ( |3.6| ) and (3/7) to this solution, 
and taking the brane tension to be T{(/)) = Tie"^*^, we find 

30/i 
and 

32 ± 16Ub - Uj - IWj ^ ^^M^^^^(x-2)<^. 9^ 

VV2^ 6(l±?7fe)V2 Vl2^rl{x-2)^/T±Ul 

A few comments are now in order. If = then x = 2 and therefore Ti must be 
fine-tuned [^]. Otherwise fixing x-, the dilaton coupling to the brane, fixes the value of 
the integration constant ^ for given A. Then, the brane tension parameter Ti in eq. ( |6.9| ) 
still depends on the arbitrary integration constant in eq. ( |6.4[ ), which could conceivably 
tune itself to allow four-dimensional flat solutions for any values of the tension. This is the 
basic idea behind 'self-tuning'. It is also easy to see that the brane tension can be positive 
for both branches. 

If we relax the assumption of Z2-symmetry across the brane, the parameter may 
take different values on each side of the brane, and so we will have one extra arbitrary 
integration constant in the expressions for x a-nd Ti. 



- 15 - 



6.3 Stability and four-dimensional effective gravity 

The functions appearing in the graviton equations (0|-46|) are: 



w 



=FV3a 



r2 50U^U^ - 



12C/V1 ± Ub 



P = T- 



(l±[/)(C/2 + 5C/2 



(6.10) 



(6.11) 



(6.12) 



where Uq = U{(j. = 0) = y^l — A/ Ac and the lower sign corresponds to the Einstein branch. 
In these coordinates r ranges from to oo. We will only consider the cases without naked 
bulk singularities {Uq > for the Gauss-Bonnet branch, and < C/q < 1 for the Einstein 
branch) . 

First note that the bulk graviton equation ( [4.10 ) is singular in the limit r — > cxd, since 



the functions w ( 6.10 ) and p ( 6.11 ) blow up. This is always true except for the trivial 



case of the constant zero mode. In the Schrodinger formulation ( 4.16 ), this translates into 
a potential V{x) diverging at a finite distance from the brane in the rc-coordinate. The 
situation is thus similar to that of a particle confined in a box, and we shall impose that 
the wave-functions vanish "at the edge of the box" . This is equivalent to require: 



~ 



] for 



oo 



(6.13) 



All the modes satisfying eq. ( |6.13| ) are normalisable on the bulk interval with respect to 
the weighting function w{r) ( |6.10| ). Considering the asymptotic behaviour of the wave 
equation ( 4.10| ) with the expressions ( |6.10 ) and ( 6.11| ) shows that eq. ( 6.13D selects a one- 
parameter family of solutions for ipm on each side of the brane. Continuity on the brane 
and the junction conditions ( [4.11 ) then result in a homogeneous system of two equations for 
the two unknown parameters. Vanishing of its determinant then gives the allowed values 
for m?, which are therefore discrete. Since a is the only other mass scale appearing in the 
problem, we expect the order of magnitude of mass spacings between the modes to be: 



Am ■ 



a 



(6.14) 



Secondly note from eq. ( 6.11[ ) that the sign of p is always negative for the Gauss-Bonnet 
branch. This implies that the sign of the kinetic term in an effective action for 7^1/ is also 
negative. This branch therefore suffers for a graviton ghost in the bulk. The same is true 
for the corresponding original six-dimensional black-hole vacuum [E^], although it has been 



argued to be classically stable |47|. On the other hand, p is always positive for the Einstein 
branch, so this case is ghost free. As discussed in section |^, there may be tachyon modes 
if either w or [v] is negative. This could occur for the Einstein branch, in which case the 
solution would be unstable. 
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To illustrate the above two points consider the case // = 0. We can then find analytic 
expressions for the modes. Indeed define: 



A 



\m\- — -^^r^e 



For > the bulk solution satisfying ( 6.13| ) is 



r 



and the mass spacings between the Kaluza-Klein states are given by: 



12a 



(6.15) 



(6.16) 



(6.17) 



in agreement with ( |6.14 ). 

Now consider a solution with m? < 0, satisfying ( |6.13 ) at infinity, 



^™ ;572 ^5/2 I — 



(6.18) 



In the case of two Einstein branches on each side of the brane, the corresponding junction 



condition ( 4.11 ) then implies 



5/2 



3C/n 



4(1 - Uo) 



3/2 



(6.19) 



If < C/q < 1) this equation is always satisfied by one real value of A_ , and so there is always 
a tachyon mode when fi = for the Einstein branch, even though p and w are positive. 
For small Uq, m? ~ —[UQ/a)e^^''/rf^. The higher order curvature term has destabilised 
the solution. Note in contrast that the junction conditions do not allow for any tachyonic 
gravitons in the case where the Gauss-Bonnet branch is present on at least one side of the 
brane (and = 0). 

To ensure that u; > we require 5(5 — 2-v/6)C/q < < 5(5 + 2\/6)C/q. To ensure 
[v] > we need Ui, < (4/3)(l — \/l — 15Uq/16), and so fi must be negative. It is only 
possible to satisfy both these conditions simultaneously if Uq > 0.87, hence A/Ac must be 
small if the spacetime is to be free from tachyons and ghosts. 

As for the constant zero mode, we see from eq. (6.10) that the volume element ( [4.15| ) is 
infinite for the solution (6.3-6^), in agreement with the general discussion of section ^. In 
the present example this is because, in order to avoid naked singularities, we have kept the 
region of space-time which contains the conformal infinity of the original six-dimensional 
space ( |6.1D , which is asymptotically anti-de Sitter. Nevertheless, it is possible to localise 
gravity on the brane, if we break Z2 symmetry by joining together an Einstein with a 
Gauss-Bonnet branch. The price to pay however is a bulk graviton ghost on the Gauss- 
Bonnet side. Indeed, for = 0, the solution ( 4.14 ) of the perturbation equation ( 4.1C| ) 



which satisfies the boundary condition ( 5.13 ) is given by 



5m 



In 



^ + - 1 ± In 
6C/2 6 ' 



l±U 



V5Uo arctan 



V5{U-Uo) 



(6.20) 
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for /i 7^ and C/r^ if // = 0, where the constant C may be different on both sides of the 
brane. The junction condition ( 4.11| ), together with continuity of V'o on the brane, then 



require 7^ and two different branches of the solution on each side of the brane. The 
special junction conditions for higher order terms discussed in section ^, as well as the two 
possible signs of p (and hence the ghost) are essential here. In this case, the zero-mode 
(6.20) is normalisable, and separated from the massive ones by the mass gap ( |6.14 ). Since 



it is natural for a to be of the order of the fundamental scale of the underlying theory, 
the contribution of the massive gravitons to brane gravity at low energies is likely to be 
negligible. This setup is still free of singularities in the bulk, and the tension of the brane 
is not fine-tuned. 

7. Conclusions 

In this paper we have investigated the properties of a Poincare brane world, located in a five- 
dimensional background containing a scalar field. We considered the most general quartic 
order action which produced field equations that were linear in second order derivatives, 
and contained no higher derivatives. If the scalar field plays the role of the dilaton, the 
couplings of the dilatonic higher order terms in the Lagrangian ( |2.lD are constrained in 
order to reproduce the scattering string amplitudes on shell |^4|. We have also derived 
their values (2.7) when the scalar field arises as the unique modulus coming from the 



Kaluza-Klein toroidal compactification of extra dimensions. 

In order to study how gravity may be perceived by a four-dimensional observer, we 
considered background perturbations corresponding to graviton modes which are tangent 
to the brane. Although the tensorial structure of bulk gravitons is still the same as in 
Einstein gravity, the higher order terms lead to significant modifications to the original 
mechanism of gravity localisation. In particular, they lead to mixed rather than Neumann 



boundary conditions ( 4.11 ), much like induced gravity models 0, where the momenta of 
modes appear in the boundary conditions themselves. This renders questions of eigenmode 
orthogonality and of background stability far more subtle and constrained than the Einstein 



case (|4.12 ). For example it is possible to have a scalar field background which is stable 



in Einstein gravity, but which develops tachyon modes when a Gauss-Bonnet correction 
is included. A solution with this type of instability ( |6.18| ) was described in the previous 
section. 

Furthermore, for some non Z2-symmetric configurations, gravity may be localised on 
the brane, even though the volume element is infinite. This involves a non-constant gravi- 
ton zero mode, and is only possible if there are multiple branches of solutions, as occurs in 
Einstein-Gauss-Bonnet gravity. However, combining different branches in this way generi- 
cally implies a bulk graviton ghost on one side of the brane. 

In section |5| we presented a method for analytically obtaining the general solution of 
the field equations, which as we showed in section ^, is a necessary step if one is to address 
the 'self-tuning' mechanism in this context (i.e. trying to resolve the four-dimensional 
cosmological constant problem by allowing four-dimensional Poincare invariance on the 
brane to be obtained without fine-tuning of its tension, due to its conformal coupling to 
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the scalar field ||30|]). We then argued that requiring finite volume element, as well as no 
naked singularities (and finite string coupling), generally implies fine-tuning rather than 
self-tuning of the brane tension. 

Finally, to illustrate all of the above, we considered in detail the case where the scalar 
field is the modulus of a sixth fiat dimension. The general solution then is simply ob- 
tained by analytic continuation of the corresponding six-dimensional Gauss-Bonnet black- 
hole The spatial topology of the black hole horizon transforms into the four- 
dimensional spacetime curvature of the brane thus providing the de Sitter, Poincare, and 
anti-de Sitter brane solutions. Regular brane solutions can then be constructed without 
fine-tuning of tension, or naked singularities in the bulk. It is rather intriguing that the 
zero-mode graviton is generically non-normalisable i.e. gravity is not localised, except if we 
consider different branches of solutions on each side of the brane. Such a setup allows for 
a normalisable zero-mode graviton to be separated from the (discrete) massive ones by a 
mass gap, of the order of the fundamental scale coupling of the higher order terms. The 
drawback of this case is the presence of a bulk graviton ghost in the Gauss-Bonnet branch 
which signals a potential quantum inconsistency of the gravitational background. This is 
a rather subtle problem since the corresponding original six-dimensional black-hole anti-de 



Sitter vacuum |4^] has also been argued to be classically stable |47] to small perturbations 
due to the change of sign of the AD energy in the Gauss-Bonnet branch. The presence, 
significance and eventual cure of such bulk ghosts certainly demands further study. 

Furthermore, on a different note, higher dimensional toroidal compactification should 
relate the scalar field solutions discussed here (for > 2) to the Gauss-Bonnet equivalent 



of cosmic p-brane solutions [ 52 1 . Such solutions yield backgrounds of higher co-dimension 
brane worlds ||53| , It would be intriguing to know their Gauss-Bonnet versions and 
to see if they can solve some of the problems of their Einstein counterparts such as the 
presence of the bulk naked singularity or the fate of self-tuning in higher co-dimension 
spacetimes (see for example ||55[|). 
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A. Appendix: Field equations 

Varying the action ( |2.lD with respect to the metric gives: 

fAgab + fGab - f'[VaVb(^ " gabDc/)] - /" [Vb(/>Va(/> - <?afe(V</))2] 

- 4a;/ (VbCf>VaCf> - ^gab{^<pf^ + 2a {fHab + 2Paebc{f"^'<P^'4> + /'V^V^0)) 

- 8ba {f[2V'<PV''<PPaebc + {^(pfGab + 2(V, Ve0) (VfeV^^) - 2iVa 
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- Va</.Vb0n0 - (V0)2VaVfc0 + 2v^0V(„(/>V6)Ve0]) 

+ \f[{^<l^?9ab - 2(V0)2v„</)Vb<; 

- 16ca/ (^2(V(/>)2Va,/.Vb0 - ^(V0)^(7a6) = (A.l) 
where the prime denotes differentiation with respect to 0, and 

Hah = RadecRh^^^ — '^^R^^Reacb — '^■RaeR^h + RRab — ^ffafe-^GB (A. 2) 

Paebc = Raebc + Rbedac + Racdbe — Recdab — RabQec + -^RgabQec — -^RgbeQac ■ (A. 3) 



Varying {2A) with respect to the scalar field gives: 

= -8a;/n0 - 4ujf'{V<pf - f'R + 2f'A + 2A' f + a (-/'£gb - 16(3c/' + af") (V0)^ 

- 166 G"'' {f'Va(bVb(p + 2fVaVb(p) + 32a/ {Dcpf - 64c/ (V^)^^^ 

- 64(2c/ + af) vVV^VfeVa^ - 32a/ (V^V^VaVb,^ + iJ^V^^Vfe^)) (A.4) 
For a conformally flat ansatz of the form (|3.l| ), the and (5,5) components of 



|A.l|) read, respectively: 

4(/." - 6A" - 4(w + 2)(j)'^ - 12A'^ - 2Ke^^ + 12 A' (j)' + 6A'B' - 4.(t>' B' 

+24ae-2^ (^'2^" - 2A'20" - 40'^'^" + 4A'20'2 - Qcj,' A"" + A'^ - ^'^5' + dA'^E'cj)') 

+486 ae~2S (^'^A" + 2^W" + 2^'^^'^ - 2A! (f)"^ - 3A'B'(I)'^) 

+32aae"2^ (-,^.'2,/," + ^'^ + B'<p''^) - IGcae-^^,/.'^ = (A.5) 
e2S2A _ 4u;(/.'2 + 12^'^ - 160'^' 

+ (24[-A'4 + SA'^^'] - 2886 [^'^^'^J + 32a [4^^'^ + 0'^] - 48c [(/)'^]) = (A.6) 



while the scalar field equation ( A.4 ) gives: 
-4w(/." - 8^" + 4w(/.'2 - 20A'2 - 2Ae2^ - lQujA!(t)' + 8A'B' + Aujcj)' B' 

+ae-2^ (24[4^'2^" + 5A'^ - AA^^B'] - 966 [A'V + 2(1)' A' A" - A'^c))'^ + 40' A'^ - 3A'2(/>'5']) 

+32a (2</.'2^" + 2</.'20" + 4A',/.V" + 8A'V^ - ^'^ - 6A'B'<P'^ - 2B'(^'^) 

-16c (6(/.'2(/)" + S^V^ - 34>'^ - 65>'^) = (A.7) 

where a prime denotes derivative with respect to z. 
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